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ABSTRACT 

Two sets of experimental data were examined 
in this paper, ground and space experiments, 
for cylindrical canisters with thermal energy 
storage applications. A 2-D computational 
model was developed for unsteady heat transfer 
(conduction and radiation) with phase-change. 
The radiation heat transfer employed a finite 
volume method. The following was found in this 
study: 1) Ground Experiments, the convection 
heat transfer is equally important to that of the 
radiation heat transfer; Radiation heat transfer 
in the liquid is found to be more significant than 
that in the void; Including the radiation heat 
transfer in the liquid resulted in lower 
temperatures (about 15 K) and increased the 
melting time (about 10 min.); Generally, most of 
the heat flow takes place in the radial direction. 
2) Space Experiments, Radiation heat transfer 
in the void is found to be more significant than 
that in the liquid (exactly the opposite to the 
Ground Experiments); Accordingly, the location 
and size of the void affects the performance 
considerably; Including the radiation heat 
transfer in the void resulted in lower 
temperatures (about 40 K). 

1 . Nomenclature 

A surface area 

c material specific heat 

D c directional cosine 

e enthalpy 

g gravitational acceleration 

H m PCM heat of fusion 

h heat transfer coefficient 

I radiative intensity 

l b blackbody intensity 

k thermal conductivity 

Nu Nusselt number 

q heat flux 

r radial coordinate 

r position vector 

S radiative source function 
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Sm 

modified radiative source 
function 

s 

direction vector 

T 

temperature 

T 0 

cooling fluid inlet 
temperature 

T f 

cooling fluid temperature 

T m 

PCM melting point 

t 

time 

u 

overall heat transfer 
coefficient 

V 

volume 

V 

velocity of the liquid PCM 

Xf 

PCM quality = e / H m 

z 

axial coordinate 

z rad 

z-coordinate of the 
radiator emitting surface 

a 

absorption coefficient 

P 

extinction coefficient 

Pm 

modified extinction 
coefficient 

P 

density 

e 

emissivity 

<D 

scattering phase function 

a 

scattering coefficient 
(Stefan-Boltzman 
constant in equation (10)) 

Q 

solid angle 

Subscripts and superscripts 

b 

bottom 

E 

finite volume to the right of 
the control volume 

eff 

effective 

L 

liquid 

m 

metal 

mz 

mushy zone 

0 

outer 

rad 

radiative 

s 

solid 

t 

top 
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2. INTRODUCTION 

Freezing and melting due to heat transfer 
take place in different engineering applications. 
At present time, phase-change phenomena 
have become an indispensable part of many 
processes, for example, manufacturing of glass 
crystals, and metal alloys, casting processes, 
food preservation and cryopreservation of 
biological cells, just to name a few. A recent 
review of the freezing and melting process was 
provided by Yao and Prusa (1 989). On the other 
hand there are several engineering applications 
which involve energy demand at a later period 
than its supply, which requires a thermal storage 
system. The choice of an appropriate storage 
system is most frequently between using a 
Sensible Heat Storage (SHS) material and 
Phase-Change Material (PCM). A review of the 
two systems has been provided by Beasley and 
Clark (1984). nvestigations of PCM systems 
have been provided by Beasley et al (1989), 
Kerslake and Ibrahim (1993), & (1994) and 
Adebiyi (1991). Recently more experimental 
data have been obtained for ground PCM 
system, Kerslake (1991), as well as Thermal 
Energy Storage (TES) experiments conducted 
on the Shuttle Columbia. These efforts were 
part of the development of Solar Dynamic 
Systems in space and were provided for 
different canister design, material, dimensions 
and type of salt. 

Several studies were conducted to utilize the 
available experimental data for both the ground 
testing and the space ones, in order to provide 
a better understanding for the phase-change 
process. Kerslake (1 991 ) simulated the ground 
experiment while Namkoong et al (1995) 
simulated the TES-1 space experiment. The 
motivation of this study is to further enhance 
the computational capabilities and physical 
understanding of the phenomena by including 
a proper modeling of the void in the canister. 
Further more, the radiation heat transfer 
analyses, in both the void and the PCM liquid, 
were improved by using a finite volume method 
Chai et al (1994). In this paper, a systematic 
approach has been taken for the modeling 
process. The predictions of the temperature 
history will be obtained with no radiation, 
radiation in the void, radiation in the void and 
the liquid. The predictions will be evaluated 
against the experimental data. One of these 
data sets were conducted at NASA Lewis 
Research Center (Ground Experiments), the 
other was part of the OAST-2 Hitchhiker payload 
on the Columbia Shuttle STS-62 in early 1994 


(Space Experiments). Detailed information 
about the Ground Experiments can be found in 
Kerslake (1991), while those for the Space 
Experiments can be found in Namkoong et al 
(1995). The physical characteristics of the 
canisters used in the two experiments are given 
in Table I. 

3. Numerical Analysis 

3.1 Geometry and Boundary Conditions 

3.1 .1 Ground Experiments 

The temperature and density distributions in 
the system are assumed to be axisymmetric 
which reduces the problem to a 2-dimensional 
cylindrical geometry (r,z). The geometry of the 
problem is shown on Figure 1 . The void is 
considered to be in the region adjacent to the 
top wall, consistent with the experimental 
observations when the canister was in the 
vertical position. The volume fraction of the void 
is 30% of the total canister volume. The void’s 
volume and shape are assumed to be constant. 

The top and the bottom walls of the canister 
are insulated. Heating is provided by a heater 
on the outside wall.The heat flux is uniform and 
periodic in time. The effective heating power 
which is the difference of the heater power and 
the power of the radiative heat loses between 
the canister and the test chamber was 
computed by Kerslake (1991) and shown in 
Figure 2. 

Cooling of the canister is provided by cooling 
air, flowing in the inside tube. The air inlet 
temperature is also provided by Kerslake 
(1991). 

3.1.2 Space Experiments (TES1) 

The geometry of the TES1 experiment is 
illustrated on Figure 3. Since the TES1 
experiment was performed in a microgravity 
environment, the location of the void is 
controlled by the surface tension rather than 
the gravity force. Because of the Marangoni 
effect, the void tends to be in the hottest region 
of the canister that is the outer lower corner in 
this problem (Namkoong et al. (1995)). The 
void’s volume fraction is 15%. 

Again the top and the bottom sides of the 
canister are insulated. Heating is provided by 
a heater on the outside wall of the canister. 
Similar to the ground experiment, the effective 
heating power can be estimated (Figure 4). 
Cooling is done by radiation from the surface 
of a circular disk (radiator) attached to the 
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canister via a conductor rod. The boundaries 
of the rod and the radiator are insulated except 
for the radiating surface of the radiator. 

In both experiments, the objective is to 
determine the temperature distribution in the 
system. The following four cases are 
considered: (1) neither the void nor the liquid 
PCM are transparent to the thermal radiation; 
(2) the void is transparent but the liquid PCM is 
not; (3) both the void and the liquid PCM are 
transparent. (4) same as case (3) but the liquid 
PCM absorbs the thermal radiation with an 
absorption coefficient of 0.84 cm' 1 . 

The similarities and differences of the 
temperature distributions in the four cases are 
of a particular interest along with evaluation of 
the computational results against the 
experimental data. 

3.2 Governing equations 

PCM, walls, conduction rod and radiator 
The enthalpy formulation of the heat transfer 
equation is used. This approach has proved to 
be effective for problems with solid-liquid phase 
change (Kerslake and Ibrahim (1 993), Wichner 
et. al. (1 988), Solomon and Wilson (1 986)). The 
conservation of energy equation in this 
approach becomes: 

(37 + v • V)(pe) = V • (kVT) - V - q nd (1) 

dr 

This approach eliminates the need to introduce 
an infinitely large heat capacitance in the 
melting region to insure the isothermal melting. 
The enthalpy is connected to the temperature 
via the constitutive equations: 

T m +e/c s : e <0 (solid PCM) 

T= T m : Otte#H m (mushy PCM) (2) 

Tm+fe-HJ/CL : H m <e (liquid PCM) 
T m +e/c m : -4#e#4 (metal parts) 

Indeed, the temperature stays constant and 
equal to the melting temperature in the so-called 
Amushy zone® which can be viewed as a 
transition zone between the solid and the liquid 
phases. The second term of equation (1) is a 
transport term and it contains the velocity of 
the liquid phase. The velocity distribution in the 
liquid PCM is governed by the gravity and 
surface tension of the liquid-void interface. 
Determining this distribution is complicated and, 
most importantly, time consuming. Since we 
are interested mostly in the global heat transfer, 
the natural convection in the liquid PCM is 
modeled using the existing empirical 
correlations. The common approach is to 


estimate the enhancement of the value of the 
liquid PCM thermal conductivity such that: 

kL,eff=Nuk L (3) 

The liquid PCM will be shown to be present 
in the problems considered in this study in the 
form of a vertical layer. Therefore, the 
correlation for a vertical layer with isothermal 
or isoflux boundary conditions by Ozisik (1 985) 
is used. 

The effective conductivities of the solid PCM 
and the metal parts are equal to their normal 
values, the liquid effective conductivity is 
computed using equation (3), and the effective 
conductivity of the mushy PCM is treated as a 
linear function of the PCM quality Xf=e / 

kmz.eff = (1 * Xf) kg + Xf ki if} ff (4) 

The liquid PCM (LiF) has been found 
(Wichner (1987)) to be transparent to the 
thermal radiation for wave-lengths <5.5 pm and 
semi-transparent for wave-lengths >5.5 pm. 
The average Rosseland absorption coefficient 
is virtually constant in the temperature interval 
of interest (1000 K - 1200 K) and equal to 
0.84 cm' 1 (Wichner (1 987)). It is of one’s interest 
to examine both the transparent and the semi- 
transparent models of the liquid PCM. Note 
that in the former case the last term of equation 
(1) becomes zero, since it can be shown that 
the divergence of the radiative heat flux is 
proportional to the absorption coefficient for a 
uniform medium (see also equation (1 4) below). 
However, the radiative heat flux can be non- 
zero in the case of the transparent medium, 
although its divergence is zero. Therefore, this 
heat flux has been computed and used in the 
energy balance equations for the solid PCM and 
walls finite volumes adjacent to the liquid PCM 
layer. 

Void 

The void is assumed to be filled with the LiF- 
CaF2 vapor in the Ground Experiments and the 
LiF vapor in the Space Experiments with 
negligible thermal capacity and absorption 
coefficient. Hence, the heat transfer equation 
in the void is 

V.(/cVT) = 0 (5) 

Similarly to the liquid PCM model, the radiative 
heat flux in the void can be non-zero, although 
its divergence is zero due to the negligible 
absorption coefficient. Therefore, this heat flux 
has been computed and used in the energy 
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balance equations of the PCM and walls finite 
volumes adjacent to the void. 


radiative intensity l(r,s) and the source function 
S(r,s) (Modest (1993)): 


The following three boundary conditions in the 
experiments have been applied: 

(a) The “heating” boundary condition at the 
outside wall: 


r ~ fo • 'km - Qin ( 6 ) 

(b) The insulation boundary condition at the 
top and the bottom walls: 

z = zt , Zb : -km = 0 ( 7 ) 

(c) The “cooling” boundary conditions are 
different in each experiment. 

(cl) Ground Experiments: 

r = n : -k m ^~ = u(T- Tf(z)) (8) 

dr 

The cooling air temperature is determined by 
utilizing the time-dependent inlet profile 
(Kerslake (1991)) and assuming a linear 
temperature distribution in the z-direction. 

(c2) Space Experiments: 

r= 0 : ^ = 0 (9) 

dr 

Z=Zrad : -km^ = <Stm(T 4 -Tj ink ) (10) 
Tsjnk is the temperature of the surrounding. 


i» miw vvi muivi i 

Initially, the system is assumed to be at a 
uniform temperature in both experiments. The 
initial temperature is 922 K in the Ground 
Experiments and 283 K in the Space 
Experiments. 


3.3 Numerical method 

The simple explicit method is used to solve 
equation (1). However, the heat transfer 
equation in the void (5) is elliptic in contrast with 
the PCM heat transfer equation, which is 
parabolic. Therefore, a simple explicit method 
cannot be used, and the successive over- 
relaxation (SOR) by lines procedure, for the 
Laplace equation (Anderson et. al. (1984)), has 
been utilized . 


= - P(r)/(r,s) + S(r,s) (11) 
os 

where p(r) = <x(r) + <j(r) is the extinction 
coefficient, a(r) is the absorption coefficient, a(r) 
is the scattering coefficient, and 

S(r, «) = a(r) l b (r) + J l(r,s'Ms',s)dCl' 

An 

( 12 ) 

In these equations, r is the position vector and 
s is the directional vector. If equation (11) 
together with emissive boundary conditions 
(see Modest (1993)) is solved then one can 
compute the radiative heat flux and its 
divergence via the following equations: 

q(r) = J l(r,s)sdQ (13) 

An 

V-g(r) = a(4 nl b (r)- J /(r,s)dft) (14) 
An 

To solve equation (11) numerically, the 
coordinate space is discretized into finite 
volumes in the usual manner. Since the intensity 
/ depends also on the angular direction, the 
angular space should be discretized as well. 
In the present method, the angular space is 
divided into solid angles. 

Integrating equation (11) over a control finite 
volume AV and a control solid angle Aft 1 and 
using the divergence theorem, we obtain 

J J i 1 (s' • n)dAdD.' = J J (-p l'+S')dV<£l' 

AO' AA AO' AV 

(15) 

where AA is the surface bounding the control 
volume AV. Following the FVM, the intensity is 
assumed to be constant within the control 
volume and the control angle. Hence equation 
(15) becomes in a 2-D axisymmetric system: 

j (s* • ni)d ft 1 = (- ^ + Sm) A I/Aft ' 

M AQ 1 

(16) 


The Finite Volume Method (FVM) developed 


by Chai et al. (1 994) has been used to compute 
the radiative heat flux. The equation of radiative 


heat transfer can be written in terms of the 


where the summation is performed over the four 
faces of the control volume, and the modified 
extinction coefficient and the source function 
are: 
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Sm = “It * -h £♦' '' 4£l1 '' 


471 ^ 

r*i 


One can solve equation (16) for the intensity in 
the center P of the control volume 

(1 - Ar/2 rp)teD l cr l' w + Arp£ r 's + s(nP ArAzAn ' 
' P (1 + A r/2 rp)AzD l cr + ArD & + p' m p ArAzACl 1 


(17) 

where the subscripts “P”, “W", and “S” denote 
the points {/,/}, {/-?,;}, and {i,j-1) respectively, and 

d' c i = J {s' nfidn 1 

An 1 

is either the r- or the z-directional cosine. 

Chai et al. (1994) modified equation (17) by 
exponentially extrapolating the intensity at the 
neighboring nodes to the sides of the control 
surface AA. Their procedure increases the 
accuracy of the scheme and eliminates the 
possibility of a physically incorrect negative 
intensity. The resulting scheme (Chai et al. 
(1994) has been used to solve for the radiative 
intensity and the heat flux. 

4. Results and Discussion. 

4.1 . Ground Experiments. 

4.1 .1 Convergence test. 

Figure 5 shows the comparison of the 
outside wall temperature predictions from the 
solutions using (10x10), (16x16), and (20x20) 
grids ((r,z) in the interior canister region). 
Significant differences in the temperature 
distribution are noted between grid sizes 
(10x10) and (16x16) while little difference is 
found between grid sizes (16x16) and (20x20). 
Grid size (16x16) has been used in this study. 

4.1 .2 Comparison with experimental data. 

Figure 6 shows the comparison of the outside 
wall temperature predicted by the numerical 
solution (four cases) with the experimental data. 
The numerical solutions capture very closely 
the important features of the temperature 
distribution, such as the slope during the one- 
phase solid heating (0-15 min), and cooling 
(65-91 min), point of the onset of melting (the 


distinct slope change at about 15 min), the 
temperature arrest during PCM melting(15-40 
min) and freezing (55-65 min). A data feature 
not well predicted by the model is the 
temperature during one-phase liquid heating 
(40-55 min). This could be attributed to a 
number of reasons: the axisymmetric 

assumption might have led to underprediction 
of the time needed to fully melt the PCM, the 
cooling fluid film coefficient and the liquid PCM 
Nusselt number might have been 
underpredicted, the void’s size and location 
were changing in the actual experiment while 
the numerical model assumed a static void. 

4.1.3 Effect of natural convection in the 
liquid PCM 

Convection usually lowers the thermal 
resistance of the fluid thus improving the heat 
transfer and lowering the temperature gradients. 
This can be seen from the comparison of the 
wall temperatures in two situations when the 
natural convection is assumed to be present or 
not present (Figure 7). The temperature in the 
absence of convection is considerably higher 
than the temperature when the natural 
convection is present, especially at the time 
interval when the liquid PCM is present (15- 
65 min). Therefore, the liquid Nusselt number 
is high during the melting stage but drops 
considerably in the freezing period. The 
negative heat fluxes at the inner and the outer 
walls during the freezing period of the cycle 
(Figure 2) cause layers of solid PCM to form in 
the regions adjacent to the inside and the 
outside walls. The liquid PCM is trapped 
between these two layers. The temperatures of 
these layers are close to the melting 
temperature and thus to each other, hence the 
Nusselt number is close to 1 . 

4.1 .4 Effect of radiation through the void 

Because of the location of the void, its 
presence and the processes inside it should 
not be a big factor in the overall heat transfer. 
This is obvious from the comparison of the 
temperature histories of the outside wall when 
the radiation through the void is present and 
not present (Figure 6, cases (1 ) and (2)). The 
radiative heat transfer does lower the 
temperature as one would expect (since the 
additional heat flux lowers the thermal 
resistance of the void), but this improvement is 
not significant (about 2-5 K) because of the 
location of the void in this situation. While the 
processes in the void do not affect the global 
heat transfer, the void’s presence can cause 
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the formation of local non-uniformities of 
temperature (“hot spots"). As seen on Figure 8a, 
distinctive cold and hot spots are formed on the 
inside and the outside walls of the canister, 
respectively, at locations where the void is 
present. 

4.1.5 Effect of radiation through the liquid 
PCM 

Since the liquid PCM forms a vertical layer 
in the canister, the additional heat flux through 
this layer should improve the overall heat 
transfer in the system thus lowering the 
temperatures in the canister. Since the layer is 
vertical, as opposed to the horizontal 
configuration of the void, the effects of the 
radiation through the liquid PCM should be 
expected to be more significant than those of 
radiation in the void. It is obvious from Figure 5 
that the maximum temperature of the outside 
wall is lowerwhen the liquid PCM is transparent 
to the thermal radiation than when it is not by 
about 18 K. Compared to the 2-5 K improve- 
ment by the void radiative heat flux, this is a 
much more significant effect. An important effect 
of the radiation through the liquid PCM and the 
reason for a lower wall temperature is a longer 
time needed to melt all the PCM in the canister 
(see Figure 5). Since the canister thermal 
resistance becomes smaller due to the radiative 
heat flux, more energy is transferred from the 
heater to the cooling fluid through the canister. 
Hence, a smaller portion of the incoming energy 
is spent on heating and melting of the PCM, 
thus a longer time is needed to fully melt the 
PCM. This portion of energy spent to melt the 
solid PCM becomes even smaller when the 
forming liquid PCM absorbs the thermal 
radiation. This results in an even longer melting 
time and lower temperatures as seen on Fig- 
ure 5. The energy absorbed and stored by the 
liquid PCM is then released during the freezing 
part of the cycle which results in a slightly longer 
freezing time and temperatures. 

Overall, the main effect of the radiation 
through both the void and the liquid PCM is a 
higher heat flux through the canister and a more 
uniform temperature distribution in the canister. 
Comparison of temperature contours in the 
canister for case (1) (Figure 8a) and case (4) 
(Figure 8b) provides a good illustration of this 
effect. Not only the overall temperature in the 
canister is lower in case (4), but also the 
temperature gradients and the temperature of 
the hot spots. 


4.2. Space Experiments. 

4.2.1 Comparison with experimental data 

The same 4 cases discussed above have 
been considered in the Space Experiments 
simulations. Figure 9 indicates that the 
temperature history has the same features as 
the one of the Ground Experiments. However, 
the thermal arrest in the Space Experiment is 
not as pronounced as in the ground experiment 
which can be seen from comparison of Fig- 
ure 5 with Figure 9. This is a result of the differ- 
ence between the two experiments in the 
boundary conditions in the heating period. 
While the power input and output were kept 
constant in the Ground Experiments (Figure 2), 
the power output in the space experiment is a 
decreasing function of time during this period. 
Hence the difference of the input power and 
the output power is a growing function of time 
for the Space Experiments. This results in a 
higher energy stored in the canister per unit 
time and therefore a steeper rise of the overall 
temperature of the canister. 

An opposite effect is observed during the 
cooling period. The thermal arrest in the Space 
Experiments is more pronounced than that in 
the Ground Experiments. This is a result of the 
power input in the Ground Experiments being 
negative in the cooling period while it is always 
positive in the Space Experiments. 

4.2.2 Effect of radiation through the void 

The presence of radiation in the void has a 
more significant effect in the Space 
Experiments than in the ground one. The reason 
for that is a larger surface area of the void is 
exposed to the incoming heat flow. The effect 
of radiation in the void on the temperature 
distribution in the canister is seen on Figure 9. 
As expected, the outside wall temperature is 
significantly lower in the presence of thermal 
radiation in the void. This improvement is about 
30-40 K compared to 2-5 K for the Ground 
Experiments. 

The void is again the location of a hot spot in 
the canister. As seen on Figure 10a, there is 
only one hot spot (compared to the hot and cold 
spots in the Ground Experiments) adjacent to 
the canister outside wall. 

4.2.3 Effect of radiation through the liquid 
PCM 

During the heating part of the cycle, the PCM 
starts to melt first in the region close to the 
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canister bottom wall because of the significant 
z-component of the heat flux and the presence 
of the void adjacent to the outside wall. 
Therefore, the liquid PCM exists in the canister 
in the form of a horizontal layer as can be 
deduced by the isotherms in Figures 10a, b. 
Computations showed that the radiation 
through the liquid PCM, whether absorbing or 
not, does not have an effect as important as 
the radiation through the void. This is obvious 
from the comparison of the temperature 
histories for cases (3) and (4) with that in case 
(2) on Figure 9. The reason for this phenom- 
enon is the horizontal configuration of the liquid 
PCM compared to the vertical one for the void. 
This results in the heat mostly transferred 
through the void rather than the liquid PCM. 

5. SUMMARY AND CONCLUSIONS 

Two sets of experimental data were examined 
in this paper: Ground Experiments and Space 
Experiments, for a cylindrical canister with 
thermal energy storage applications. A 2-D 
computational model was developed for 
unsteady heat transfer with phase change. An 
important feature of the model is using a Finite 
Volume Method (FVM) for the radiation heat 
transfer. Three different grid sizes, 10X10, 
16X16 and 20X20 have been attempted. The 
later two cases showed that the numerical 
solution is grid independent. The predictions of 
the temperature history were obtained for four 
different cases regarding the radiation heat 
transfer: 1 ) no radiation in the void or the liquid 
PCM, 2) radiation in the void only, 3) radiation 
in the void and the liquid PCM simultaneously 
and 4) same as case (3) except that the liquid 
PCM has a thermal radiation absorption 
coefficient of 0.84 cm' 1 . The predictions were 
evaluated against the experimental data. The 
following is summary of finding in this study: 

1. Ground Experiments: 

1 .a The convection heat transfer is equally 
important to that of the radiation heat 
transfer; l.b Radiation heat transfer in 
the liquid is found to be more significant 
than that in the void; l.c Including the 
radiation heat transfer in the liquid resulted 
in lower temperatures (about 15 K) and 
increased the melting time (about 10 min.); 
I.d Generally, most of the heat flow takes 
place in the radial direction; l.e The hot 
spots were found (as expected) in the outer 
region of the void. 


2. Space Experiments: 

2.a Radiation heat transfer in the void is 
found to be more significant than that in 
the liquid (exactly the opposite to the 
ground experiments); 2.b Accordingly, the 
location and size of the void affects the 
performance considerably; 2.c Including the 
radiation heat transfer in the void resulted 
in lower temperatures (about 40 K); 2.d 
The heat flow takes place in both the radial 
and axial directions; 2.e The hot spots 
were found (as expected) in the outer region 
of the void. 
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Table I. Test Canister Design Features, Ground Experiment and Space Experiment 


Feature 

Ground Experiments 

Space Experiments 

Material of Construction 

Haynes Alloy 1 88 

Haynes Alloy 188 

Outer Diameter/Wall Thickness 

4.98/0.129 (cm) 

6.96/0.102 (cm) 

Inner DiameterVWall Thickness*7Rod Diameter 

2.07/0.261 (cm)/N/A 

3.81/0.102/3.59 (cm) 

Length / Side Wall Thickness 

2.43/0.091 (cm) 

6.94/0.102 (cm) 

Phase Change Material (PCM) 

LiF-20 %CaF 2 

LiF 

PCM Melting Point / Heat of Fusion 

1042 K/ 816 J/g 

1121 K/ 1037 J/g 

PCM Mass / Total Canister Mass 

53/ 137 (g) 

265 / 543 (g) 

Cycle Time (% time Heating) 

91.1 min (60%) 

131 min (58%) 


* - Cooling air tube ** - Thickness of canister inner wall plus cooling air tube wall 
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Figure 1. Geometry of the Ground 
Experiments with the Void Model. 



time of the cycle, min 

Figure 4. Power Input and Output for the 
Space Experiments. 
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Figure 2. Power Input and Output for the 
Ground Experiments. 
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Figure 5. Computational Canister 
Temperature Data with Different 
Grids, Ground Experiments 
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Figure 3. Geometry of the Space Experiments 
with the Void Model. 
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Figure 6. Experimental (symbols) and 
Computational Canister Temperature Data, 
Ground Experiments. 
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Figure 7. Effect of Natural Convection on 
Canister Temperature, Ground Experiments, 
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Figure 8a. Computational Canister 
Temperature Contours, Ground Experiments, 


time = 50 min, Case (1). 
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Figure 8b. Computational Canister 
Temperature Contours, Ground Experiments, 
time = 50 min, Case (4). 


Figure 9. Experimental (Symbols) and 
Computational Canister Temperature Data, 
Space Experiments. 



r, cm 


Figure 10a. Computational Canister 
Temperature Contours, Space Experiments, 
time = 76 min, Case (1). 
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Figure 10b. Computational Canister 
Temperature Contours, Space Experiments, 
time = 76 min, Case (4). 
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